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ABSTRACT 

We present a simple method for measuring cosmological bulk flows from large redshift sur- 
veys, based on the apparent dimming or brightening of galaxies due to their peculiar motion. It 
is aimed at estimating bulk flows of cosmological volumes containing large numbers of galaxies. 
Constraints on the bulk flow arc obtained by minimizing systematic variations in galaxy luminosi- 
ties with respect to a reference luminosity function measured from the whole survey. This method 
offers two advantages over more popular bulk flow estimators: it is independent of error-prone 
distance indicators and of the poorly-known galaxy bias. We apply the method to the 2MASS 
redshift survey (2MRS) to measure the local bulk flows of spherical shells centered on the Milky 
Way (MW). The result is consistent with that obtained by Nusser and Davis (2011) using the 
SFI++ catalogue of Tully-Fishcr distance indicators. We also make an assessment of the ability 
of the method to constrain bulk flows at larger redshifts (z = 0.1 — 0.5) from next generation 
datasets. As a case study we consider the planned EUCLID survey. Using this method we will be 
able to measure a bulk motion of ~ 200 kms -1 of 10 6 galaxies with photometric redshifts, at the 
3cr level for both z ~ 0.15 and z ~ 0.5. Thus the method will allow us to put strong constraints 
on dark energy models as well as alternative theories for structure formation. 



Subject headings: Cosmology 



Introduction 



In the standard cosmological paradigm, devi- 
ations from the Hubble flow, commonly named 
peculiar motions, are the result of the gravita- 
tional force field dominated by the dark matter, 
with luminous galaxies behaving like test particles. 
Hence observations of the peculiar velocity field of 
galaxies are a direct probe of the three dimensional 
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distribution of the dark matter. Cosmological bulk 
flows of spherical volumes centered on the MW are 
the most intuitive moments of this velocity field. 
These bulk flows are typically measured over large 
spheres (i.e. £ 30ft. -1 Mpc) where linear grav i- 
tational instability holds (e.g. iNusser et ah 1991 ). 
facilitating comparisons with cosmological models. 

So far, estimates of the bulk flow in the lo- 
cal Universe within ~ 100/i -1 Mpc from us have 
mostly relied on galaxy peculiar velocities inferred 
from known distance indicators and measured red- 
shifts. The distance indicators commonly used to 
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estimate peculiar velocities arc based on well de- 
fined, if often heuristic, relations between intrin- 
sic, observable galaxy properties, one of which de- 
pends on the galaxy distance . The typical example 
is the Tully-Fisher relation (jTullv fe Fisherlll977l ) 
between the absolute magnitudes and the rota- 
tional velocity of spiral galaxies. Distance indi- 
cators have been extremely useful in enforcing our 
confidence in gravitational instability as the main 



mechanism for structure formation ( Davis et al 



201l[ ) and in putting constraint s on cosmologi- 
cal mo dels. A recent analysis bv iNusser fc Davisl 
(|2011ah has shown that the bulk flow inferred 
from a trimmed version of the SFIH — h catalog 
of spir al galaxies with I-band Tully-Fisher dis- 
tances (iMasters et all 120061: ISpringob et al.l 12007 : 
Davis et al.l 2011 ) is consistent with the standard 
ACDM cosmological model with the best fit pa- 
rameters of the Seven- Year Wilkinson Microwave 



Anisotropy Probe (WM AP7) (e.g. iJarosik et al 



20101: lLarson et al.ll2010l). Furthermore, using the 
recovered bulk flow, INusser fc Davisl (|2011al) de- 
rived an interesting constraint on the growth rate 
of fluctuations in the Universe at z ~ 0. 

The paucity of distance indicators available and 
the observational difficulty in measuring the rele- 
vant quantities make peculiar velocity quite diffi- 
cult to estimate. In addition, their accuracy de- 
grades linearly with redshift, limiting their useful- 
ness to rather small rcdshifts (cz < 10000 kms^ 1 ). 
More accurate peculi ar velocity measurements 
based on Tyla SN fe.g. lRiess et al.lll997l : lDai et al 



2011) and surface brightness (e.g. iBlakeslee et al 
19991 ) are possible for substantially fewer galaxies. 
As a result, galaxy peculiar velocities have been 
estimated for a relatively small number of galax- 
ies. To make the matter worse these galaxies are 
not uniformly distributed across the sky, poten- 
tially inducing unwanted biases in the bulk flow 
estimate. Here we propose to circumvent these 
potential problems induced by the use of distance 
indicators and constrain the bulk flow using only 
galaxy luminosities and redshifts. 

Rcdshifts of galaxies systematically differ from 
the actual distances by the line-of-sight compo- 
nent of the bulk motion. Hence, the intrinsic 
luminosities of galaxies inferred from the observed 
flux using redshifts rather than distances redshifts 
appears to be brighter or dimmer. In presence of 
large bulk flows the effect is systematic and its 



strength depends on on amplitude and direction 
of the bulk motion. Object-by-object magnitude 
variations as a result of peculiar motions are over- 
whelmed by the natural spread in the distribution 
of magnitudes. Hence, individual peculiar veloci- 
ties cannot be derived from these considerations. 
However, an estimate of the bulk motion of some 
sobvolume within the survey can be derived on a 
statistical basis, by comparing the luminosity dis- 
tribution of galaxies in the subvolume with the one 
in the whole survey. 



Th e method has a long history. iTammann et al 



(Il979h correlated the magnitudes with rcdshifts of 
galaxies to constrain the velocity of the Virgo clus- 
ter relative to the Local Group (LG) of galaxies. 
The main limitation of this method has been the 
limited number of the objects available and the 
limited size of the sampled volume. The extension 
we propose here is timely in view of the consid- 
erable progress in current and future redshift 
surveys. 

In §2 we introduce the method, presenting gen- 
eral expressions and deriving the relevant approx- 
imations. In §3 we apply the method to the Two 
Mass redshift Survey (2MRS) of - 23000 galaxies 
limited to magnitude K = 11.25. In §4 we discuss 
prospects for successful applications of the method 
to future redshift surveys. We conclude in §5 with 
a general discussion. 

2. The method 
2.1. The set-up 

Consider a subvolume in a survey of galaxies 
with measured redshifts cz (in kms -1 ) and ap- 
parent magnitudes m limited to m < mi. The 
redshifts are given in the frame of reference of the 
cosmic microwave background (CMB) radiation. 
We envisage two descriptions for the geometry of 
the subvolume. (i) a local subvolume of a thick 
shell centered on the observer (the MW) in an all 
sky redshift survey, and(ii) a distant region where 
all galaxies closely lie along the same line-of-sight. 

Let r (also in kms -1 ) be the (unobserved) lu- 
minosity distance to a galaxy. For simplicity of 
notation and description we assume here that the 
distance and spatial extent of the survey are small 
so that r is well approximate by the physical dis- 
tance. The results can readily be extended to the 
general case once we specify the underlying cosmo- 
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logical model. The unknown, true absolute mag- 
nitude, M, of a galaxy is expressed as 



acquires the simpler form 



M 



15 — 5 log r = Mq — 7 



(1) 



where the measurable absolute magnitude Mo = 
m — 15 — 5 log cz is determined from observations, 
and 7 = 51og(r/cz). The observed redshift of the 
object is cz = r + v where v is the line-of-sight 
component of its peculiar velocity. A bulk mo- 
tion of the subvolume yields a systematic differ- 
ence between r and cz which induces a mismatch 
between M and Mo . Therefore, the bulk motion of 
the subvolume (relative to the motion of the whole 
survey) can be estimated by demanding that the 
distribution of measured magnitudes, Mq, is con- 
sistent with the distribution of true magnitudes, 
M, in the whole survey. This is the underlying 
principle of the method outlined here. 

For the case (i) of a local thick shell centered 
on the MW, a bulk flow B yields a galaxy radial 
peculiar velocity v B = B cos 6 where 6 is the an- 
gle between B and the line-of-sight to the galaxy. 
This flow introduces a systematic angular dipolar 
modulation in the distribution of Mq — M across 
the sky, which allows a determination of the mag- 
nitude and direction of B. For the case (ii) of a 
distant region, only the component of B in the di- 
rection to the line-of-sight to the subvolume is rele- 
vant, giving rise to a systematic difference Mq — M 
for all galaxies in the subvolume. In this case 
only the line-of-sight component of B can be con- 
strained. 

Systematic differences between Mq and M can 
be appreciated by comparing the luminosity distri- 
bution of galaxies in the subvolume, i.e. their lu- 
minosity function, with the luminosity function of 
the whole survey. We define the luminosity func- 
tion, $(M), expressed in terms of the absolute 
magnitudes, as the number density of galaxies per 
unit magnitude. We assume that $(M) is well ap - 
proximated by a Schechter form ( Schechter 1980f ) 



$(M) = 0.41n(10)$*10°- 4(Q+1)(M *- M) 

X CX P (_l °-4(M*-M)^ (2) 

The normalization $* and the shape parameters 
M* and a generally depend on the galaxies' type, 
redshift and band of observation. In terms of the 
luminosity (M — — 2.5 log L + const), this function 



= 0.41n(10)$* ^j-^j exp [^J^j ■ 

* (3) 

2.2. The General Formalism 

We quantify the effect of the bulk flow via the 
probability distribution functions of the magni- 
tudes in the subvolume and the whole survey. We 
write the conditional probability that a galaxy at 
redshift cz in the subvolume possesses a measured 
absolute magnitude Mq, 

P(M \cz;v B ) = J P{M Q \r)P(r\cz)dr 

= J P(Mo\M)P(M\M int )P{M int )P{r\cz) 

x6(Mi (r) - M)dM int dM dr . (4) 

Here v B is the component of the bulk flow, B, in 
the line-of-sight to the galaxy. The magnitudes 
Mq and M are related by Eq. Q]and M, is the in- 
trinsic absolute magnitude which differs from M 
due to photometry errors and small scale peculiar 
motions not described by the bulk flow. The Heav- 
iside step function, 0, accounts for the magnitude 
cut imposed by the apparent magnitude limit mi, 
i.e. M;(r) = mi — 15 — 51ogr. The intrinsic 
magnitude, M lriti , appears only in the underlying 
luminosity function P(M- ln t\r) and in P(M\M- m ti) 
which accounts for the difference between M and 
Mi n t. The probability distribution function 

P(M)= f P(M\M int )P(M int )dM int (5) 



is proportional to the luminosity function of the 
whole survey, $(M), and we assume that it can 
be described well by a Schechter form despite the 
convolution of P(M int ) with P(M\M int ). All other 
terms in Q are straightforward. We have: 



P(M \M) = 8 U {M + 7 - M ) 



(6) 



where S D is the Dirac 5-function. P(r\cz) can be 
written in a more convenient form using Bayes' 
theorem: P(r\cz) — P(cz\r)P(r)/ P(cz) where 

P(r) = r 2 n(r) and P{cz) = (cz) 2 n(cz) , (7) 

and we assume that the number density of objects 
is constant along the line of sight, n(r) sa n(cz) ~ 
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const , and where 



P(cz\r) 



(8) 



which assumes that redshifts are normally dis- 
tributed about the value r + v B . The dispersion is 
a = <Jq + a 2 where 00 corresponds to the rms of 
errors in the measured redshifts and a v to small 
scale motions not described by the bulk flow. Sub- 
stituting all this in Eq. |4] and integrating over M 
gives 

/>oo 

P(M \cz;v B ) oc / r 2 dr$(M -7) 



x Q(Moi(cz)-M )e 2 ^ (9) 

where the argument of the step function is now 
M i(cz) - M where M m {cz) = M/(r) + 7 = m - 
15 — 5 log cz. 

Instead of integrating over radial coordinates, it 
is convenient to integrate over luminosity L(r) oc 
10 -o.4Af(r)_ Substituting r = cz(L/L ) 1 ' 2 and us- 
ing the Schechter functional form ([3]) for $(M(L)) 
we get 

/>oo 

P{M \cz;v B ) ex Lq 3/2 / dLL 3 / 2+a e~ L ^ L * 

Jo 

The conditional probability must integrate to 
unity, hence we write the normalized probability 
for the observed Mq as 



P(M a \cz;v B ) 



0.41n(10) Ll +a J °° dye F ^ 



where 

F(y)= (3 + 2 a )lny-2/ 2 L /L 



(11) 

(y + i; B /cz-l) 2 



2((T /cz) 2 

(12) 

and y = (L/Lq) 1 / 2 . Note that the expression does 
not involve the G function. This is because the 
measured Mq is derived from observed redshifts 
and apparent magnitudes and hence it is guaran- 
teed that 6 (M[(cz) - M ) = 1. Equation [TT] is 
the expression for the distribution of the observed 



magnitudes of galaxies in the subvolumc. By com- 
parison, the expression for the distribution of the 
observed magnitudes in the whole survey is 

P SU rvey(M)|cz) = J P(M \cz; V B )P(v B )dv B . 

(13) 

For a Gaussian field P(v B ) is expected to be Gaus- 
sian with zero mean and rms o~b ■ The integration 
over v B gives a similar expression to (|lip but with 
v„ = and a 2 = a 2 + a\ instead of a 2 . 

D C z cz D cz 

The general strategy to estimate a bulk mo- 
tion should now be clear. It can be described 
as a two-step porcedure. (i) Find the parame- 
ters a and of the Schechter function which 
minimize the quantity — ^ - lnP survey ((Moj \czj), 
where the summation is over all galaxies in the 
parent survey, (ii) Insert these parameters in () 1 1 1) 
and search for the value of v B which minimize 
— l n -P(M)i|czi; w B ), where now the summation 
is only over galaxies in the subvolume. 

The terms (a cz /cz) 2 and (a cz /cz) 2 are very 
small compared to the expected signal |B|/cz, 
even in presence of photometric errors, so the ap- 
proximation is valid for all applications considered 
in this work. For |v B |/cz,<C 1 and a^jcz <C 1 
the function F(y) is very well approximated by a 
quadratic function in y. This fac t allows us to ap- 
ply th e steepest descent method (jArfken fc Webed 
120051 ) to evaluate numerator and denominator of 

(ED. 

For the numerator we write F(y) = F max + 
(l/2)(y - y m ax) 2 /Cr 2 , where y max is t he value 
where F(y) has a maximum and Cj= = 1/ \f— d 2 F/dy 2 
is also evaluated at y = y m ax- The steepest de- 
scent method yields the following approximation 



/>oo 

/ dyc F = V^Cjre 1 
Jo 



(14) 



For the denominator we integrate over Lq and 



get 



/ dL L% / dye F = L\ +a / dyt 

JLi(cz) JO JO 



G(y) 



where 



G = lnT{l + a,y 2 Li/L*)+lny 



(15) 

(y + V B /CZ- 1) 2 



^cz/czf 



(16) 
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where T(a, x) = J°° dte~ t t a ~ 1 is the upper incom- 
plete gamma function. The application of steepest 
descent to the denominator of Eq. [15] gives the fol- 
lowing expression for the conditional probability: 

P(M \cz;v B ) =0Aln(W)(j^j + ^e F — " G — 

(17) 

where G max is the maximum of G and Cg = 



d 2 G/dy 2 is evaluated at this maximum. 
These approximations are very accurate. We 
have checked that for the typical values which are 
relevant for this work the accuracy is better than 
one part in 10 . 

2.3. Long distance and small redshift er- 
rors approximation 

As a final step we now assume that the am- 
plitude of the bulk flow, small-scale peculiar ve- 
locities and redshift errors are small compared to 
to the typical redshift of the objects in the sur- 
vey. In this hypothesis we can use (fT7j) to derive 
P(A1q\cz;v b ) to first order in v B /cz and a cz /cz. 
The maximum of F, to second order in <j rz jcz 
and v B /cz, is obtained at y max = I — v B /cz + 
(3 + 2a — 2L n /L r )(cr cz /cz) 2 . Substituting this ex- 
pression into (|12[) we find, to first order terms in 
v/cz, that 

e F max = f i _ ( 3 + 2a )^ e -(i-2«B/c»)W£. , 

(18) 

and similarly 

e G — = (l - ^) T (1 + a, (1 - 2v B /cz)Lt/L*) . 

(19) 

Further, we get Cg 2 /Cjr 2 = 1. Inserting all these 
expression into (JTTJ) gives 



P{M Q \cz;v B ) = 



0.41n(10) (fc) 



1+a 



r(l + a,Zi/i a 



(20) 

where Lq = (1 — 2v b /cz)Lq and Li = (1 — 
2v B /cz)Li. In practice, the effect of a bulk flow 
vb is that of shifting the estimated luminosity of 
a galaxy Lq by a factor which is proportional to 
the amplitude of the flow. If applied to a distant 
subvolume moving at velocity vb, the net effect 
would be that of shifting the typical luminosity 



L* of the measured luminosity function, with no 
impact on the slope a and normalization $». This 
result is hardly surprising, but it is reassuring to 
see it obtained in the limit of small velocity \v/cz\ 
and errors a cz /cz. Second order corrections to 
this expression involve terms of 0((a cz /cz) 2 ) and 
0((v B /cz)) 2 , which are typically very small when 
one considers cosmological volumes. 

3. Application to 2MRS 

In this section we apply the method outlined 
above to the 2MRS. The goal is to estimate the 
bulk flow of a thick spherical shells centered on 
the MW at z ~ 0. 2MRS is an all-sky redshift 
catalog of about 23,200 galaxies, which is com- 
plete down to the magnitude K = 11.25. De- 
tails about the catalog, including the precise com- 
pleteness, sky coverage and sele ction effects can be 
found in (jHuchra et al.ll2005bl ). The preparation 
of the catalog for the purpose of t he application of 
the me thod is is done similarly to Nusser fc Davisl 
(|2011ah . 

The description of the method presented in 
the previous Section assumes that the luminosi- 
ties of galaxies in the sample are drawn from a 
single luminosity function. However, the method 
could easily be extended to account for the differ- 
ent luminosity functions of different galaxy types 
(e.g. red vs. blue, spirals vs. ellipticals) and 
environments (high density vs. low density re- 
gions). For this reason we have divided galaxies in 
two samples: spirals and ellipticals, and measured 
their luminosity functions separately. We find 
(a,M») = (-0.803,-23.53) for early-type and 
(a,M») = (-0.888,-23.12), for late-type galax- 
ies in agree ment with t he pr evious, independent 
estimates by IWestoverl (|2007T ) . We made no at- 
tempt to determine $* which is not relevant for 
our goal. These are the reference values for the 
luminosity function of the whole survey. To de- 
termine the bulk flow we consider a spherical sub- 
volume of radius cz = 10000 kins -1 that contains 
16460 galaxies brighter than K = 11.25, of which 
10366 and 6094 are late and early types, respec- 
tively. We divided this sample into spherical shells 
each one 4000 kms -1 thick. To estimate the three 
Cartesian components of the bulk flow, B, in each 
shell, we have minimized the probability function 
— Si ^P(M 0i \cZi] v B ) with respect to B, where 
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the summation is over all galaxies in the shell. The 
expression for (|20|) is used for P where different 
values of a and L* are used for late and early type 
galaxies. The results are shown in Fig. [T] The 
points represent the Cartesian components of the 
bulk flow, indicated on the plot, estimated in cor- 
respondence of the mean radius of each shell. The 
la errorbars were obtained from 200 mock sam- 
ples mimicking the 2MRS subvolumc. Each sam- 
ple contains a number of points similar to that of 
the 2MRS galaxies. The points are randomly 
distributed within a sphere of 200/i~ 1 Mpc and 
are assigned absolute magnitudes according to the 
overall reference luminosity function of all 2MRS 
galaxies. Mock galaxies were assigned random 
peculiar velocities sampled from a Gaussian dis- 
tribution with zero mean and a 300 kms -1 scatter 
to mimic the combined effect of small scale veloc- 
ity dispersion and errors in the measured redshift. 
No bulk flow was assigned to the mock samples. 
Our minimization procedure was then applied to 
the 200 mock samples to derive the intrinsic un- 
certainties in the estimate of the bulk flow. These 
uncertainties are plotted as la errorbars. 

Despite the large errors on the bulk from 
K=11.25 2MRS, the results are very encour- 
aging. They provide a strong incentive for a 
more thorough study using the deeper K=12.2 
2MRS which will triple number of galaxies within 
cz = 10000 kms -1 . The results can b e com - 
pared with bulk flow of Nusser &; Davis! ( 2011a ) 



(see their figure 5) obtained using the Tully-Fisher 
SFI++ data. A full comparison between the pho- 
tometric and the SFI++ results is complicated 
and should take into account the complex covari- 
ance of errors in both measurements. The effort 
is redundant at this stage due to the large er- 
rors on the current 2MRS photometric bulk flow. 
Therefore, we only make a comparison between 
the two bulk flows at R « 60/i _1 Mpc. At this 
intermediate radius the photometric bulk B p hot ~ 
(100±90, -240±90, 0±90) kms" 1 and the SFI++ 
B S fi++ ~ (15 ± 20, -280 ±30, 120 ± 20) (se e top 
panel of figure 5 in iNusser fc Davis] ( 2011ah ) are 
fully consistent. We note that B p hot refers to 
shells of thickness 40/i _1 Mpc, while BgFi++ to 
the full sphere within R. But this hardly matters 
since B p hot versus distance is nearly flat within 
the errors. The reason that B p hot is computed 
for shells rather than the whole spheres is that 



the bulk is heavily weighted by the lower redshifts 
galaxies due to the rapid decrease in the number 
of 2MRS galaxies at larger redshfits. Comput- 
ing the photometric bulk in spherical shells al- 
lows a clearer inspection of the constraints we get 
at larger redshifts. In contrast, the SFI++ flow 
within a sphere of radius R is, however, mainly 
dominated by galaxies nearer to R. 

The method outlined above is not the only one 
to estimate the bulk velocity of a spherical vol- 
ume centered on the MW. We have implemented 
an alternative procedure based on the fact that 
a bulk flow, B, introduces a dipole-like angular 
modulation in Mq across a spherical shell cen- 
tered on the observer. For |B|/cz <C 1, Mo = 
M — 51oge B ■ r/cz = — 2.1715B • r/cz, where f 
is the usual unit vector in the linc-of-sight to the 
galaxy. The bulk flow, B, can then be found by 
minimizing x 2 = XXMh + 2.1715B • r,) 2 . We 
have checked that this alternative strategy gives 
results that are noisier than but consistent with 
the luminosity-based method which uses the shape 
of the luminosity function, in addition to the mea- 
sured magnitudes and redshifts. 

4. Prospects for future surveys 

In this section we assess the potential of our 
method applied to future datasets. The main 
goal of most of the next-generation large red- 
shift survey is to constrain Dark Energy proper- 
ties. For this reason, they are designed to op- 
timize observations at z ~ 1, i.e the epoch at 
which Dark Energy started dominating the energy 
budget of the Universe. However, some of these 
projects will also provide redshifts, measured ei- 
ther spectroscopically or through multi-band pho- 
tometry, for a large number of relatively local 
galaxies (i.e. z < 0.5), where the bulk flow can 
be reliably estimated by our photometry-based 
method. Among the planned redshift surveys that 
meet th ese constraints are the ground-based LSST 
project ( LSST Science Collaborations et al.j 2009h 
and the EUCLID satellite mission (|Laureiidl20QgF . 
The two surveys enjoy the same large sky coverage 
(~ 20000 deg 2 ) and are expected to deliver pho- 
tometric redshifts for a large (2-4 billion) number 
of galaxies in the redshift range z = [0,3]. In ad- 
dition, the EUCLID project is also expected to 
measure spectroscopic redshifts for about 70 mil- 
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Fig. 1. — The bulk flow of 40ft- 1 Mpc-thick shells centered on the Milky Way. The X, Y and Z components 
of the bulk flow arc plotted against the mean radius of the shell, as indicated in the figure. 



lion galaxies. Given the present uncertainties in 
the survey design we do not attempt here to pro- 
vide an accurate forecast for the possibility of mea- 
suring our bulk flows from these surveys. Instead, 
we will allow for generous variations in the survey 
design parameters and provide an order of magni- 
tude estimate for our capability to constraint the 
bulk flow at z < 0.5. For this purpose we will use 
the EUCLID survey as a baseline and allow for 
significant variations in the number of the objects 
in the survey which should include the case of the 
LSST survey. 

To model the characteristic of both redshift sur- 
veys, the spectroscopic and the photometric ones, 
we need to specify the luminosity function of the 
galaxies, their redshift distribution and the errors 
on the photometry and the redshift. 

For the luminosity function one we consider a 
Schechter form in the H AB-band with va l ues o f 
a and A/* derived from Cirasuolo et al. ( 2010l ) 



assuming the H-K col o r corr ection proposed by 
J lmencz fc Kashlinskvl f)l999h . The luminosity 
function evolves with redshift. Here we are in- 
terested in two epochs: z = 0.1 and = 0.5, for 
which we find (a,Af„) = (-1.07,-22.15) and 
(a,M„) = (-1.07, -22. 21), respectively. The spec- 
troscopic catalog is expected to contain objects 
with apparent magnitude brighter tha Hab = 19.5 
while the photometric catalog will sample galaxies 
as faint as Hab = 24. We assume a photometric 
error a m = 0.2 which represents the expected er- 



ror for faint objects and thus overestimate the un- 
certainties for brighter objects. We note that the 
expected number of objects in the EUCLID photo- 
z catalog (~ 10 5 gals, deg -2 with Hab < 24) 
is remarkably similar to that in the LSST one 
(~ 2 x 10 5 gals, deg" 2 with i AB < 25.3). A 
similarity that also applies to the expected galaxy 
redshift distribution and that further justifies our 
choice of not considering the two surveys sepa- 
rately. 

The galaxy redshift distribution dN/ dz for the 
E UCLID spectroscopi c catalog has been modeled 
bv lCeach et ail f)2010h . In the current design, the 
survey will target Ha galaxies above a limiting 
flux /ho = 3x lO" 16 erg cm^s" 1 . The dN/dz 
of these objects is listed in Table 2, column 3 of 
that paper. The numbers must be further divided 
by a factor ~ 4 to account for the efficiency in the 
determining the redshift from the Ha line (Guzzo, 
private communication). This survey will measure 
the spectroscopic redshift of objects with z > 0.45 
with an expected accuracy a z — 0.001. In the 
redshift bin Az = 0.1 centered at z = 0.5, the 
expected number of galaxies is N ~ 4 x 10 6 . To 
simplify our calculation we adopt a more conser- 
vative N = 10 6 and also explore the pessimistic 
case of measuring only N = 10 5 redshifts. 

The expected number of galaxies in the photo- 
z catalog can be derived by the observed galaxy 
dN/dz the GOODS and UKIDDS-UDS fields 
(Zamorani, private communication). This cata- 



log will include low-redshift objects and redshift 
will be measured with an accuracy cr z /(l + z) = 
0.05 — 0.03 (expected-goal). Here we adopt the 
first, more conservative estimate. At z = 0.5 
the number of objects in this catalog is ~ 10 
times larger than in the spectroscopic one, i.e. 
the same shell at z = 0.5 considered above will 
have N = 10 7 objects with measured photo-z. 
As before, we will also consider two pessimistic 
cases in which the number of measured rcdshifts 
is N = 10 6 and N = 10 5 . In the following we will 
be also interested in measuring the bulk flow in 
the nearby universe, at z = 0.13. At this distance 
EUCLID is expected to measure N = 10 6 photo-z 
in a sphere of radius 100 h^ 1 Mpc. We consider 
this as a reference number density at z ~ 0.1 but 
also consider a pessimistic and an optimistic cases 
in which N = 10 5 and N = 10 7 , respectively. 

To investigate whether our method can be ap- 
plied successfully to these next-generation surveys 
we consider two possible estimates of the bulk flow: 
at z = 0.13 and at z = 0.5. As a first example we 
consider a bulk flow v B = 300 kms -1 in a sphere 
of radius 100 h^ 1 Mpc at z = 0.13, similar in size 
to the one we have considered in the 2MRS sam- 
ple. As a second case we consider a large bulk 
flow vb = 1000 kms" 1 of a shell at z = 0.5 with 
Az = 0. 1, reminiscent of t h e rece nt claimed detec- 
tion by iKashlinskv et al.l (|201dh and potentially 
conflicting with the standard cosmological frame- 
work. These two bulk flows introduce the same 
systematic shift in the apparent magnitude of the 
objects: 6M w 0.018. Our task is to assess if the 
next generation surveys described above will be 
able to detect such shift. To summarize, we will 
explore the following scenarios: 

• The case of a spherical region of radius 100 
h~ x Mpc moving at v B = 300 kms -1 at 
z = 0.13, containing N = 10 6 objects with 
measured photo-z. In addition to this refer- 
ence case we also consider an optimistic and 
pessimistic scenarios in which the number of 
measured photometric redshifts is N = 10 7 
and N = 10 5 , respectively. The results rela- 
tive to this scenario are listed in the first 12 
rows of Table [TJ 

• The case of a shell at z — 0.5 centered 
on the MW, with Az = 0.1 and moving 
v B = 1000 kms" 1 , sampled by N = 10 7 ob- 



jects with measured photo-z (reference case). 
We also consider the pessimistic cases of 
N = 10 6 and N = 10 5 objects. The results 
relative to this scenario are listed in rows 13- 
24 of Table Q] 

• The case of the same shell at z = 0.5 
moving at the same speed but sampled by 
N = 10 6 with redshift measured spectro- 
scopically. We also consider a more pes- 
simistic case of iV = 10 5 objects. The results 
are listed in rows 25-26 of Table [1] 

The general strategy to assess our ability in 
measuring the bulk flow is straightforward. First, 
we consider a very large set of objects represent- 
ing the whole survey to which we assign abso- 
lute magnitudes according to the reference lumi- 
nosity function and determine the a and M* by 
minimizing — P survey (Mo J \cz). Then, for each 
of the three scenarios itemized above, we repeat 
the exercise considering 70 subsamplcs, i.e. we 
determine a set of (a, A/*) values by minimizing 
— J2j P(Moj\cz; vb = 0). Objects in the subsam- 
ples are assigned a magnitude error by adding a 
Gaussian scatter of zero mean and width gm and 
no bulk flow. The condition for a significant bulk 
flow detection is that the rms scatter in the value 
of M* measured from the 70 mock samples is larger 
than the expected magnitude shift 6M w 0.018 

The results of this test are summarized in Ta- 
ble [T] They refer to the three scenarios outlined 
above. To assess the impact of photometric errors 
we have explored the case of ctm = (no pho- 
tometric errors) and that of fixing a at its mean 
value as derived from the mean of subsamples. In 
this Table, [i SM is the difference between the mean 
of the recovered M* from the 70 subsamples and 
the true value and its scatter a SAU . The quantities 
fi^ and <7 fo arc defined similarly but for a. 

For N = 10 7 , we get very small a SM and <7 fa . 
This means that one can treat the corresponding 
values as those that would be obtained from the 
full survey (which is expected to contain a signif- 
icantly large number of galaxies than 10 7 ). Since 
the presence of a bulk flow is not expected to af- 
fect a but only M*, this result suggests that one 
can safely fix a to the value measured from the 
whole survey, i.e. that the case which is directly 
relevant for our method is the one tagged 'fixed a 
in Table [TJ We also show the case of free a for the 
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sake of completeness. The results also show that 
when (tm = the values of n SM and vanish. 
This is a known effect and reflects the fact that 
random photometric errors preferentially scatter 
galaxies to the brighter end of the luminosity func- 
tion. However, the bias would affect the estimate 
of the true luminosity function from the whole sur- 
vey and hence it is irrelevant to our method which 
only uses the values of a and M* measured from 
the full survey in order to get the bulk flow by 
minimization of — P(Moj \cz; vb)- 

The results show that the bulk flows that we 
are considering here and which would introduce a 
systematic magnitude shift SM = 0.018 would be 
detected at high significance in all three reference 
cases considered (boldface quantities in Table [TJ. 
In particular, we see that the 300 kms" 1 bulk flow 
at z — 0.13 would be detected at ~ ha significance 
by measuring the photometric redshift of 10 6 ob- 
jects in a a sphere of 100 h^ 1 Mpc. At z = 0.5, 
10 7 photometric redshifts would allow us to detect 
a bulk flow of 1000 kms -1 at ~ 20a, decreasing 
to ~ 7 and 3a, respectively, for 10 6 and 10 5 photo- 
metric redshifts. The statistical significance of the 
bulk flow detection can be further increased by re- 
ducing redshift errors, i.e. by using spectroscopic 
rather than photometric redshifts. The result is 
seen by comparing the two cases of a spectroscopic 
and a photo-z with the same value of N and om, 
(e.g. row 25 vs. row 20 in Tabled]). The value 
of au is reduced by a factor of ~ 2. The reason 
of this improvement is that galaxies in the spec- 
troscopic catalog sample the brighter end of the 
luminosity function which is sensitive to the value 
of M*. The result is that, for the same value of N, 
the galaxies in the spectroscopic sample are able 
to constrain M* considerably better than galaxies 
in the photometric catalog. The same reasoning 
explains why a SM in row 16, which refers to a sam- 
ple of distant, bright galaxies at z = 0.5 is smaller 
than the value in row 4 which refers to a sample 
of nearer and fainter galaxies. 

5. Discussion 

Over the years cosmological bulk flows has 
been estimated from the measured peculiar ve- 
locities of a large variety of object s r anging from 
galax i es dGiovanelli et al.l 1998alfbl; Dekel et al 



Sarkar et al.ll2007) clusters o f galaxies (Lauer fc Postman 
1994 iBranchini et al.l Il99l Hudson et al.l |2004| ) 



and SNIa (jRiess et al.l 119951) . Conflicting results 



triggered by the use of error-prone distance indi- 
cators have fueled a long lasting controversy on 
the amplitude and convergence of the bulk flow 
that is still on For example usin g the SFI++ 
galaxy catalog ( Watkins et al. 20091 ). claimed the 
detection of a bulk flow of 407 ± 81 km s _1 within 
R = 50 h~ 1 Mpc, inconsistent with expectation 
from the ACDM model. This result has been 
challenged by the re-ana lysis of the same data by 
Nusser fe Davisl (|2011bl ) who found a bulk flow 



amplitude consistent with ACDM expectations. 
Several bulk flow estimates have been recently per- 
formed from the dipole-like anisotropy induced by 
the kinetic Sunyaev-Zel'dovich decrement in the 
WMAP temperature map, measured at the posi- 
tion of X-ray galaxy clusters. When interpreted 
as a coherent motion, this signal would indicate 
a gigantic bulk flow of 1028 ± 265 km s _1 within 
R = 528 ft. _1 Mpc, clearly inconsistent with the 
standard picture of gravitational instability. How- 
ever, a more r ecent and independe nt analysis of 



20101 ) did not con- 



WMAP data (|Osborne et al 
firm this result. Finally, systematic anisotropies 
in the observed clustering of objects induced by 
peculiar velocities, the so called redshift space 
distortions, are a popular and reliable way to esti- 
mate the growth rate of density fluctuations and 
to constr ain the amplitude of l arge-scale coherent 



motions dPeacock et al. 2001 ; Guzzo et al. 2008 ; 



Blake et al.ll2010h . However, their interpretation 



19991 ; ICourteau et al.l l200d: Ida Costa et al.l |2000|; 



in terms of bulk velocity is not straightforward. 

Here we have presented and implemented a 
method to estimate cosmological bulk flows that 
only depends on photometry and redshift mea- 
surements and not on distance indicators. The 
main requirements for a successful application of 
this statistical method is dense sampling over large 
regions of the universe. A large number of object 
is mandatory since, for a given bulk flow ampli- 
tude, errors are dominated by Poisson noise, as 
clearly indicated by the fact that the error on M* , 
a SMt , listed in Table Q] scales as V^V- We find that 
the accuracy in the measurement of galaxy magni- 
tude and redshifts are less crucial in the bulk flow 
estimate. This fact indicates that photometric- 
redshift catalogs are better suited for this anal- 
ysis than the spectroscopic ones which typically 
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Table 1: Expected error on from a survey similar to EUCLID. Col. 2: redshift. Col 3: number of galaxies 
in the subvolume. Col 4: assumed error on the magnitude. Col 5: deviation from the actual M* written as 
the systematic bias ± the la error. For "fixed a" , the value of a was fixed from the whole survey. Rows 1-24 
and 25-26 refer, respectively, to the photometric and spectroscopic redshifts catalogs. The differences results 
between these catalogue for the same N are due to the different apparent magnitude limits: Hab = 24 for the 
photometric and Hab = 19.5 for the spectroscopic catalog. Most relevant cases for EUCLID are emphasized 
in boldface. 



row $ 


z 


N 


cm 




6a ^ Set 


1 


0.13 


10 v 





0± 0.0013 


0± 0.0002 


2 


0.13 


10 7 





± 0.0009 


fixed a 


3 


0.13 


10 7 


0.2 


-0.035 ±0.0016 


-0.0028 ±0.0002 


4 


0.13 


10 7 


0.2 


-0.035 ±0.0013 


fixed a 


5 


0.13 


10 e 





0± 0.0041 


0± 0.0007 


6 


0.13 


10 6 





0± 0.0031 


fixed a 


7 


0.13 


10 6 


0.2 


-0.035 ±0.0056 


-0.0028 ±0.0007 


8 


0.13 


10 6 


0.2 


0.035 ±0.0037 


fixed a 


9 


0.13 


10 5 





0± 0.015 


0± 0.0025 


10 


0.13 


10 5 





0± 0.0011 


fixed a 


11 


0.13 


10 5 


0.2 


-0.035 ±0.0158 


-0.0028 ±0.0021 


12 


0.13 


10 5 


0.2 


-0.035 ±0.0114 


fixed a 


13 


0.5 


10 v 





0± 0.0014 


± 0.0005 


14 


0.5 


10 7 





0± 0.00075 


fixed a 


15 


0.5 


10 7 


0.2 


-0.05 ±0.0015 


-0.01 ± 0.0005 


16 


0.5 


10 7 


0.2 


-0.05 ±0.0008 


fixed a 


17 


0.5 


10 B 





± 0.0048 


0± 0.0018 


18 


0.5 


10 6 





0± 0.0024 


fixed a 


19 


0.5 


10 6 


0.2 


-0.05 ±0.0047 


-0.01 ±0.0016 


20 


0.5 


10 6 


0.2 


-0.05 ± 0.0026 


fixed a 


21 


0.5 


10 5 





0± 0.013 


± 0.0049 


22 


0.5 


10 5 





0± 0.008 


fixed a 


23 


0.5 


10 5 


0.2 


-0.05 ±0.0136 


-0.01 ±0.005 


24 


0.5 


10 5 


0.2 


-0.05 ±0.0065 


fixed a 


25 


0.5 


10 6 


0.2 


0.08 ±0.0015 


fixed a 


26 


0.5 


10 5 


0.2 


-0.08 ±0.0044 


fixed a 
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include significantly less galaxies. Of course, this 
assumes that systematic errors and the fraction 
of 3cr outliers (which can be significant at z < 



0.5 ILSST Science Collaborations et al.l (j2009t) ) 



the photo-z measurement can be conveniently re- 
duced. Dedicated analysis that accurately mimic 
the expected performance of the surveys will be 
required to address this issue. 

Large volumes and dense sampling are require- 
ments that will be met by next-generation red- 
shift surveys like EUCLID and LSST. Therefore, 
we have carried out a feasibility study for such 
surveys. Taking into account that the sensitivity 
of our method decreases with the redshift we have 
only considered galaxies with z < 0.5. We focused 
on two different measurements: a bulk flow of 
~ 300 kms -1 at z ~ 0.15 on scales ~ 100ft -1 Mpc 
and a much larger bulk flow of ~ 1000 kms -1 
at z ~ 0.5 over a Az = 0.1 shell centered on the 
MW. Both flows would introduce the same shift in 
the measured magnitude of the galaxies. The re- 
sults, summarized in Table[T]are very encouraging. 
They suggests, for example, that one could obtain 
independent estimates of cosmological bulk flows 
(with typical amplitude of ~ 300 kms -1 ) over the 
~ 50 2MRS-sizcd samples within z = 0.2 It is 
worth stressing that photo-z are accurate enough 
for our purposes. A comparison between rows 8 
and 25 in the Table shows that 10 6 galaxies (with 
photo-z) the error in determining M* is actually 
smaller at z ~ 0.5 than z ~ 0.13. This implies that 
a bulk flow of 200 — 300 kms -1 (which amounts 
to SM ~ 0.0045 at z ~ 0.5 and S M ~ 0.012 at 
z ~ 0.13) is constrained at the 3<r level at both 
rcdshifts. 

As a preliminary application of the method, we 
have resorted to the 2MRS all sky survey which 
is limited to K = 11.25. The re sults are consis- 
tent wi th the bulk flow derived bv lNusser fc Davisl 
(|2011al ) using the SFI++ Tully-Fisher catalog of 
distance indicators. The 2MRS photometric bulk 
flow derived here is noisy due to the relatively 
small number of galaxies. Better results will be 
obtained using the deeper K = 12.2 2MRS. This 
already planned extension of the original survey 
will increase the number of observed redshifts by a 
factor of 3 within a distance of 100ft -1 Mpc. Since 
Poisson noise dominates the error budget, we ex- 
pect a factor of y/3 reduction in the size of the 
crrorbars. 



There are strong indications that the lumi- 
nosity function, $, depends on the large scale 
environment of galaxies (e.g. Balogh et al. 200 it 



Hutsi et al.l 200 



2005 



Park et al 



^ .Croton et aL. „ 
20071: iMo et alj|2004h . These studies explore the 
dependence on the density of galaxies on scales of 
a few Mpcs. As outlined above, we are interested 
in significantly larger scales of 100ft -1 Mpc. 
When averaged over these scales, the rms density 
fluctuations is less than 0.07 (for standard ACDM 
model with cr 8 = 0.8). So far no study has ad- 
dressed environmental dependences of galaxy lu- 
minosities over such large scales. If the overden- 
sity of the large scale environment is the main 
relevant parameter, then an extrapolation of the 
observed dependence on a few Mpc scales to 100 
ft -1 Mpc yields an exceedingly small effect. More- 
over, in next generation surveys we could easily 
account for such effects thanks to our ability to 
measure the bulk flow over independent volumes 
that could be classified according to their average 
density. To illustrate the robustness of our 2MRS 
result over environmental effects, we note that a 
300 kins -1 bulk flow at 60ft -1 Mpc gives a mag- 
nitude difference of SM = 0.2 between the apex 
and the anti-apex of the bulk motion. In ACDM, 
the rms density in spheres of 60ft -1 Mpc is 0.13. 
Assuming again that the overdensity of the envi- 
ronment is the only relevant parameter, we find 
that SM should dominate the environmental de- 
pendence extrapolated from smaller scales. A fur- 
ther evidence that environmental effects did not 
affect our bulk flow estimates from 2MRS is that 
the metho d here gives results tha t are consistent 
with those Nusser k, Davisl ( 2011a ) who derive the 
bulk flow based on the Tully-Fisher distant indi- 
cators, rather than photometry. 
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